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' Abstract 

> : 
b : 

. . . , Neutrino oscillation scenarios involving large angle Un — > u e oscillations 

are disfavoured in the parameter range Am 2 /eV 2 ~ by recent results 

from the CHOOZ reactor-based V e disappearance experiment. For this reason 
we extend our previous work on up-down asymmetries for various oscillation 
scenarios by computing up-down asymmetries and the R ratio for the entire 
conceivable range 1(T 4 - 10" 1 eV 2 of Am 2 . Matter effects in the Earth play 
a crucial role. We perform a \ 2 fit to the data. We find that, because of the 
matter effect, the three-flavour maximal mixing model provides a reasonable 
fit to SuperKamiokande and CHOOZ data provided that the relevant Am 2 is 
in the range 4 x 10~ 4 ~ Am 2 /eV 2 ~ 1.5 x 1(T 3 . 
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Recent data from atmospheric neutrino experiments [0 and especially the Su- 
perKamiokande experiment provide very strong evidence for large angle neutrino os- 
cillations. Traditionally the atmospheric neutrino anomaly has been represented by the 
quantity R, where 

data /-. \ 

= (N,/N e )\ MC - { > 

The quantities N e ^ are the numbers of e-like and /z-like events. In addition to an anomalous 
value for R, the Kamiokande and SuperKamiokande experiments have observed anomalous 
zenith angle dependence [||||. This zenith angle dependence can be represented by the 
up-down asymmetry parameters 



^ ~ (NaVN^)\ MC { { ) 

Here N~ v denotes the number of a-like events produced in the detector with zenith angle 
cos0 < — rj, while denotes the analogous quantity for cos 9 > r], where rj is defined to 
be positive (note that cosG > for downward going leptons). SuperKamiokande divides the 
(—1, +1) interval in cos into five equal bins. The central bin straddles both the upper and 
lower hemispheres, and is thus not useful for up-down asymmetry analyses. We therefore 
choose rj = 0.2 in order to utilise all the data in the other four bins. 

When comparing the measured results for R and the 7's with predictions from a specific 
neutrino oscillation model, the numerators are replaced by calculated predictions from the 
models, while the denominators remain as the no-oscillation predictions. Note that system- 
atic uncertainties for up-down asymmetries are expected to be smaller than for R, because 
the latter depends on the relative flux of to v e . 

The utility of using up-down asymmetries as a diagnostic tool was emphasised in Ref . [[5] , 
where up-down asymmetries were computed for various neutrino oscillation solutions to the 
atmospheric neutrino anomaly. The analysis of Ref. || focussed on the energy dependence 
of up-down asymmetries. By contrast, in Refs. |||7| we followed (super) Kamiokande and 
considered particular cuts in energy. 

In Ref. 0, we analysed four representative cases: 

(A) Maximal — v T mixing ||. 

(B) Maximal — v e mixing [ID|. 



(C) Threefold maximal mixing |TT],|T2[ amongst z/ e , and v T . 

(D) Massless neutrinos with violation of the Equivalence Principle or breakdown of Lorentz 
I3 |. The case of maximal v e — oscillations JR]] was considered for definiteness. 



mvanance 



In Ref. 0, we focussed on the region of parameter space Am 2 /eV 2 > 2 x 10 3 . However, 
large angle — > v e oscillations are now disfavoured in this parameter range because of recent 



results from the CHOOZ reactor-based V e disappearance experiment This experiment 
disfavours maximal v e — mixing for Am 2 /eV 2 > 0.9 x 10~ 3 at 90% C.L. Thus it is 
important to discuss the parameter space region Am 2 / eV 2 < 2 x 10~ 3 . In Ref we 
neglected matter effects [IIB] due to neutrino oscillations through the Earth. However, it 
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turns out that this is not a valid approximation for the 2-flavour cases for multi-GeV (sub- 
GeV) neutrinos unless Am 2 /eV 2 ~ 10~ 2 (10~ 3 ), and is never a good approximation for 
model C. Thus the purpose of this paper is to reexamine the up-down asymmetries and 
R for the entire conceivable range of interest for Am 2 (i.e. Am 2 /eV 2 > 10~ 4 ). We will 
numerically integrate the Schrodinger equation for neutrino evolution including the matter 
effects, taking the density profile of the Earth from Ref. |T7]]. We will show that for the two 
flavour —>■ v e oscillation models B and D the matter effects suppress the oscillations and 
do not improve the fit of these models to the data. Interestingly, however, for the 3-flavour 
model C the matter effects actually improve the fit of the model to the data. Although this 
model does not fit the data as well as model A, we will show that this model does provide 
an acceptable fit to the SuperKamiokande and CHOOZ data for a range of Am 2 . 

Our methodology is similar to our previous paper @ except that we have used the 
inclusive cross section given in Ref. |18| for the multi-GeV analysis. Although this cross 
section is not completely satisfactory for calculating absolute event rates because it does 
not incorporate low Q 2 effects such as the A resonance production, it is a good enough 
approximation for calculating ratios of event rates such as up-down asymmetries and R. 
We also include results for case A, even though it does not involve v e , because it will 
be interesting to compare cases A and C. [For a comparative analysis of case A and the 
somewhat similar large angle —>■ v s solution (where v s is a sterile neutrino), see Ref, 
Tj| ]. We are now able to improve on the analysis of Ref. |7[] in another respect, because we 



have been fortunate to obtain the detection efficiency functions from the SuperKamiokande 
collaboration. 

Our results are given in Figures 1-6 |2(J|2l]], together with preliminary results from Su- 
perKamiokande (corresponding to 414 days of live running) |2j 

R(sub - GeV) = 0.61 ± 0.03 ± 0.05, R(multi - GeV) = 0.67 ± 0.06 ± 0.08, 
GeV) = 0.78 ± 0.06, Y^ 2 (mulU - GeV) = 0.49 ± 0.06, 
GeV) = 1.13 ± 0.08, Yf 2 [multi - GeV) = 0.83 ± 0.13. 



Y°- 2 (sub 
Y?- 2 (sub 
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Note that only statistical errors are given for the up-down asymmetries since they should 
be much larger than possible systematic errors at the moment. 

Figures 1,2 show that all of the models A,B,C,D can provide an acceptable fit to R. 
However the up-down asymmetries Y e ^ clearly distinguish the models. The only cases which 
can provide an acceptable fit to all of the data are A and C. Indeed, the Y e ^ and R values 
for model C are quite similar to model A for low Am 2 . To understand this point, consider 
the Schrodinger equation for neutrino evolution in model C including matter effects, 



.d_ 

dx 





» " 


Am 2 








2E 
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where x is the distance travelled, E the neutrino energy, Am 2 the larger of the two squared- 
mass differences in model C and v e ^ >T {x) the wave-functions of the neutrinos. The quantity 
A(x) is related to the effective potential difference generated through the matter effect: 



A(x) 



IE 
Am 2 



V2G F N e 



2.9 x 10" 



" E/GeV " 
Am 2 /eV 2 



4 g/cm 3 



(5) 
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where Gf is the Fermi constant, N e (x) the number density of electrons along the path of the 
neutrino and p(x) the mass density of the earth (in deriving the right hand side of Eq.(|5|) 
we have assumed that the average number of protons per nucleon is approximately constant 
~ 0.48). For antineutrinos the sign of A(x) is reversed. 

For E/Am 2 values leading to large matter effects (A ^> 1/3), u e 
pressed. The system in this case exhibits approximate two flavour 
lations with 



oscillations are sup- 
— ► v T maximal oscil- 



1 — sin 



2 (L/km)(Am 2 /eV 2 ) 

3 ' (E/GeV) 



(6) 



This qualitatively explains why R, Y for case C are similar to case A for Am 2 / eV 2 ~ 3 x 10~ 3 . 
From Eq.(P), observe that for a given Am 2 the oscillation length of — > v T is not the same 
as the oscillation length for genuine two flavour oscillations - it is 1.5 times longer. This 
explains why the multi-GeV R and Y^ for model C are displaced relative to model A. Finally 
note that in the large Am 2 / E limit, the Y^ asymmetry does not approach 1 in model C. 
In fact, one can show by explicit computation that in the large Am 2 /E limit, the muon 
neutrino survival probability is related to the vacuum survival probability by 



P{v,^v tx ){A) = P{v tx ^v,){A 




1Am 2 AV 

cos — — 

3 2E 
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Note that Am 2 A/E is independent of Am 2 and E and thus it turns out that there is no range 
of parameters where matter effects can be neglected for atmospheric neutrinos in model C. 

We now perform a x 2 analysis to determine the preferred region of Am 2 for model C. We 
will not consider models B and D because they obviously lead to bad fits. We first define a 



X function for atmospheric data: [22 
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where the sum is over the sub-GeV and multi-GeV cases, the measured SuperKamiokande 
values and errors are denoted by the superscript "SK" and the theoretical predictions for 
the quantities are labelled by "th". The r\ = 0.2 choice is understood for the up-down 
asymmetries. There are 6 pieces of data in x 2 an d 1 adjustable parameter, Am 2 , leaving 5 



degrees of freedom. [Note that in the present paper we consider \U, 



a-3 I 



1/3 [a = e,/i,r, j 



1, 2, 3) for model C and thus the mixing angles do not constitute free parameters]. 

The solid line in Fig. 7 displays Xatm as a function of Am 2 for model C (also shown is 
Xatm f° r y asymmetries only). The CHOOZ experiment disfavours Am 2 /eV 



2 > 



10" 



for large angle — > v e oscillations [[15]. In order to incorporate the CHOOZ results, we 
define another x 2 '- 



Xchooz 



E 



,1 j 



5xj 



(9) 



The sum is over 12 energy bins of data (in the above equation Xj are experimental values 
from Figure 5b of Ref. [HJ and i/i are the corresponding theoretical predictions) p3[] . The 
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short-dashed line in Fig. 7 displays Xchooz as a function of Am 2 for model C. Figure 8 plots 
Xatm + Xchooz- The 3a allowed range for the atmospheric data plus CHOOZ is 



4 x 1(T 4 ~ Am 2 /eV 2 ~ 1.5 x 10~ d . (10) 

The best fit point at Am 2 ~ 8 x 10~ 4 eV 2 gives Xmin — 23 for 17 degrees of freedom 
which implies an allowed CL. of 16% for model C to explain the atmospheric data while 
simultaneously being consistent with CHOOZ. 

We have also performed a x 2 fit f° r model C to the solar neutrino data for the parameter 
space of Ref. 0]. This model predicts an energy independent solar electron neutrino flux 
reduction of 5/9 (this leads to predictions for the solar neutrino experiments which are 
reduced by a factor 5/9, except for (super) Kamiokande where neutral current effects must be 
incorporated). We have used the most recent results for Homestake (2.55±0.14±0.14 SNU) 
|H, GALLEX(76.4 ± 6.3±fJ SNU) |25[, SAGE (69.9±|;| SNU) g§. For Kamiokande 
(SuperKamiokande) we have used the 8(16) energy bins given in Ref. |S7| (Ref. 153). We 
choose to leave the boron flux as a free parameter flOjI^l and neglect the small theoretical 



errors of the other fluxes. We consider two solar models for definiteness, BP95 and TCL 
3Cfl . We find that xLiar( m ^ n ) = 44 for BP95 solar model and xLiar(™ n ) = 33 for TCL 
solar model. With the boron flux as a free-parameter, there are 27 — 1 = 26 degrees of 
freedom. Thus, the overall fit of model C to the solar+atmospheric+CHOOZ experiments 
turns out to be reasonable with Xmin = 67 (BP95) and 56 (TCL) for 43 degrees of freedom 
which corresponds to an allowed C.L. of 1% (BP95) and 10% (TCL). 

In conclusion, we have extended the analysis of Ref. to include low values of Am 2 /eV 2 . 
This is important for models which have large angle —>■ v e oscillations because of the recent 
CHOOZ results. Out of the three cases (B, C and D) which involve z/ e , only the three-flavour 
maximal mixing model (case C) provides a reasonable fit to the up-down asymmetries and 
R ratios while being consistent with CHOOZ. We thus reach the important conclusion that 
the most favoured solutions to the atmospheric neutrino anomaly in the light of CHOOZ 
are: (i) large angle or maximal — > v T oscillations (case A), (ii) large angle or maximal 



^ -> v s 



oscillations |19| and, (iii) three-flavour maximal mixing (case C) with Am 2 



8 x 10~ 4 eV 2 (sufficiently small departures from three-flavour maximal mixing would also of 
course be allowed). If (i) is true then the planned Japanese long baseline neutrino oscillation 
experiment may or may not see a signal (see Ref. [[HI), if (ii) is true then present data suggest 



that they should see a signal |TJ| , and if three-flavour maximal mixing is correct then the 
Japanese long baseline experiment has no chance to see a signal. Finally note that if model 
C is correct then the Kamland experiment should see a positive signal. 
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Figure Captions 



Figure 1. The sub-GeV R as a function of Am 2 /eV 2 [(8v /2)kmGeV] for models A (solid 
line), B (dashed-dotted line), C (dashed line) [D (dotted line)]. The usual SuperKamiokande 
sub-GeV momentum cuts have been employed. The horizontal long-dashed lines are the 
preliminary SuperKamiokande data ±lcr statistical errors. 

Figure 2. Same as Figure 1 except for SuperKamiokande multi-GeV case. 

Figure 3. The sub-GeV up-down e-type asymmetry F e 02 as a function of Am 2 /eV 2 
[(Sv/2)kmGeV] for models A (solid line), B (dashed-dotted line), C (dashed line) [D (dot- 
ted line)]. The usual SuperKamiokande sub-GeV momentum cuts have been employed. 
The horizontal long-dashed lines are the preliminary SuperKamiokande data ±1<t statistical 
errors. 

Figure 4. Same as Figure 3 except for muons instead of electrons. 

Figure 5. Same as Figure 3 except for multi-GeV sample. 

Figure 6. Same as Figure 3 except for muons and for multi-GeV sample. 

Figure 7. x 2 fit as a function of Am 2 to the SuperKamiokande atmospheric data for Model 
C. The solid line includes both R and up-down asymmetries whereas the dashed line includes 
only the up-down asymmetries. The x 2 f° r the CHOOZ reactor data (dotted line) is also 
shown. 

Figure 8. Combined x 2 fit as a function of Am 2 to the SuperKamiokande atmospheric data 
and the CHOOZ reactor data for model C. 
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